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Abstract : 

In  this  paper,  certain  path-independent  integrals,  of  relevance  in  the 
presence  of  cracks,  in  elastic  and  inelastic  solids  are  considered.  The  hyp¬ 
othesized  material  constitutive  properties  Include:  (I)  finite  and  infinit¬ 
esimal  elasticity,  (ii)  rate-independent  incremental  flow  theory  of  elasto- 
plasticitv,  and  (iii)  rate- sens  it ive  behaviour  Including  elasto-viscoplasticity , 
and  creep.  In  each  case,  finite  do  format i ons  are  considered,  along  with  the 
effects  of  body  forces,  material  .uvoler.it  Ion,  and  arbitrary  traction/displace¬ 
ment  conditions  on  the  crack-face.  Ms.  the  piivsic.il  interpretations  of  each 
of  the  inte  grals  either  in  terms  of  .  ruk-np  energy  release  rates  or  simply 
energy-rate  differences  in  two  comparison  cracked-bodies  are  explored.  Several 
differences  between  the  results  in  the  present  work  and  those  currently  con¬ 
sidered  well  established  in  literature  are  pointed  out  and  discussed. 

A 
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Introduction ; 


Since  the  pioneering  work  of  Eshelby  [1],  and  independent  discovery  of 
Rice  [2],  innumerable  number  of  papers  have  appeared  in  literature  concerning 
path-independent  integrals  (for  the  most  part  however,  concerning  the  so- 
called  J-integral)  and  their  application  to  mechanics  of  fracture.  No  attempt 
is  made  herein  at  a  comprehensive  survey  of  this  still  burgeoning  literature. 

Of  particular  relevance  to  the  present  work  are  the  important  studies,  on  con¬ 
servation  laws  in  finite  and  infinitesimal  elasticity,  by  Knowles  and  Sternberg 
[3],  and  the  interpretation  of  these  in  the  context  of  the  mechanics  of  cracks 
and  notches  in  2-dimensional  bodies  by  Budiansky  and  Rice  [4] .  It  is  noted  that 
the  studies  in  [1-4]  are  restricted  to  elasto-statics,  and  the  crack-extension 
considered,  if  any,  is  of  quasi-static  nature  i.e.,  inertia  is  considered  neg¬ 
ligible.  As  is  often  noted  in  literature,  the  so-called  J-integral  in  elastos- 
tatic  crack  mechanics  is  in  fact  one  of  the  components  of  a  vector-integral, 
and  its  relevance  is  in  the  context  of  incipient,  self-similar,  crack-extension. 

Eventhou^h  the  so-called  J-integral  [1,2]  was  intended  to  be  applicable 
to  finite  or  infinitesimal  elasticity,  its  use  had  been  extended,  in  several 
W(  rks  in  the  past  twelve  years,  far  beyond  the  range  of  its  apparent  theo¬ 
retical  validity.  Invoking  "deformation  theory  of  plasticity"  and/or  "prop¬ 
ortional  loading",  it  was  used  in  the  context  of  initiation  of  "Mode  I"  growth 
of  a  crack  in  an  elasto-plastic  body.  A  wide  class  of  literature  has  also 
grown  around  the  concept  of  the  so-called  "J^  test".  Rice  [5],  in  an  article 
which  appeared  in  1976,  and  which  is  perhaps  accurately  representative  even 
today,  succinctly  summarized  a  wide  body  of  literature  pertaining  to  the  use 
of  J  in  the  context  of  elasto-plast icity .  Later,  based  more  or  less  on 
empirical  reasoning,  J  and  the  rate  of  change  of  J  with  crack  length  ("dJ/da") 
were  postulated  to  be  "valid"  parameters  characterizing  even  stable  crack 
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growth  in  ductile  materials  [6],  wherein  crack-growth  necessarily  implies  un¬ 
loading.  However,  the  fundamental  concept  of  a  pat  Independent  integral  and 
its  physical  meaning,  if  any,  in  the  context  of  an  (even  rate-independent)  in¬ 
cremental  flow  theory  of  plasticity  yet  remains  to  be  explored. 

In  the  class  of  problems  generally  characterized  as  belonging  to  the 
domain  of  "dynamic  fracture  mechanics",  integral  relations  quantifying  the  rate 
of  energy-release  to  a  propagating  crack-tip  in  a  plane  linear  elastic  body 
undergoing  infinitesimal  deformation  were  presented  by  Freund  [7,8],  who  also 
succinctly  summarized  the  pertinent  work  of  Atkinson,  Eshelby,  Achenbach,  Sih, 
and  others.  The  path-independency,  if  any,  of  such  integrals  for  energy-release 
rates  even  in  linear  elastodynamic  crack  propagation  yet  remain  to  be  understood. 

To  the  author's  knowledge,  no  work  has  been  reported  concerning  path- 
independent  integrals,  which  may  characterize  the  severity  of  the  conditions  near 
the  crack-tip,  in  materials  cha  acterized  by  rate-sensitive  inelastic  con¬ 
stitutive  laws,  such  as,  for  instance,  visooplasticity  and  creep.  However,  in 
the  case  of  pure  steady-state  reep  characterized  by  a  power  law  (of  the  type 
e~on),  an  integral  j  (or  c*)  \  lich  is  entirely  similar  to  the  J[l,2]  for  pure 
power-law  hardening  materials  (e~on)  was  introduced  by  Goldman  and  Hutchinson  [9] 
and  Landes  and  Begley  [10),  based  on  the  observed  similarity  of  the  constitutive 
laws  in  the  respective  cases  (i.e.,  c  instead  of  e,  etc.).  However,  the 
physical  interpretation,  if  any,  of  c*  appears  not  to  have  been  fully  explored. 

The  present  work  represents  a  modest  effort  at  a  re-examination  of  path- 
independent  integrals,  and  their  relevance  to  mechanics  of  cracks,  in  elastic 
as  well  as  inelastic  solids.  The  postulated  material  behaviour  includes  the 
cases  of  finite  elasticity,  rate- independent  incremental  flow  theory  of  plas¬ 


ticity,  and  rate-sensitive  behaviour  such  as  visco-plasticity  and  creep.  In 


each  case  finite  deformations  are  considered,  along  with  the  effects  of  body 
forces,  material  acceleration,  and  arbitrary  traction/displacement  con¬ 
ditions  on  the  crack-face. 

We  start  by  considering  the  case  of  finite  elasticity  and  attempt  to 
generalize  the  conservation  laws  given  by  Knowles  and  Sternberg  [3]  to  the  case 
when  body  forces,  inertia,  and  arbitrary  crack-face  conditions  are  accounted  for.  From 
these  conservation  laws,  we  derive  a  path-independent  vector  integral  of  rele¬ 
vance  to  fracture-mechanics.  In  this  process,  we  re-examine  the  fundamental 
notion  of  "path-independence",  and  the  attendant  mathematical  and  physical 
reasoning,  as  originally  propounded  by  Rice  [2],  and  the  results  presented 
later  by  Budiansky  and  Rice  [4],  We  specialize  the  obtained  results  to  the 
case  of  linear  elastodynamic  crack-propagation.  To  understand  the  physical 
meaning  of  the  "path-independent"  integral  vector,  we  make  an  l  dependent  study 
of  the  expression  for  the  rate  of  energy-release  in  elastodynamic  crack-propa¬ 
gation.  Several  differences  between  the  results  in  the  present  v ork  and  those 
currently  considered  well-established  in  literature  are  noted  ar d  discussec  . 

In  the  second  part  of  the  paper,  we  consider  conservatior  laws,  and  the 
attendant  path-independent  integrals,  in  the  incremental  theory  of  finite- 
deformation,  rate-independent,  classical  elastoplasticity .  Once  again  we 
include  body  forces,  inertia,  and  general  crack-face  conditions  in  the  discussion. 
Also,  we  explore  the  physical  meaning  of  the  path-independent  integral  increm¬ 
ental-vector,  in  the  case  of  elasto-plasticity . 

In  the  final  part  of  the  paper  we  consider  finite  strain  rate-sensitive 
inelasticity  characterized  by  a  elasto-viscoplast ic  constitutive  law  of  the 
type  proposed  by  Perzyna  111].  We  also  treat  the  well-known  Norton's  power- 
law  tvpe  of  stcadv  state  creep  as  a  special  case.  We  point  out  certain  "incre¬ 
mental"  integrals,  which  are:  path-independent  in  a  limited  sense  in  the  case 
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elasto-viscoplastic  strains,  and  strictly  path-independent  in  the  presence 
of  combined  elastic  and  creep  strain  rates.  We  explore  the  physical  meaning 
of  these  integrals  as  well. 

Notation : 

For  convenience  to  the  reader,  we  summarize  the  notations  employed  in 
the  ptesent  as  follows:  / 


(-) 

(-) 

a  =  \.b_ 
C  =  A.  B 


A:B 


V 

~P 


Au 


t 

s 

T 

O 


under  symbol  denotes  a  vector 


under  symbol  denotes  a  second-order  tensor 

implies  a  =  A  bJ 
•  J 

implies  =  A^'B*^ 

•  K. 


implies  the  trace:  A^B 


ij 


a  second-order  tensor  in  dyadic  notation 


gradient  operator 


displacement  vector  and  its  increment 
first  Piola-Kirchhof f  stress 
second  Piola-Kirchhof F  stress 
Cauchy  stress 
Kirchhof f  stress 


a 


co-rotational  rate  of  Kirchhoff  stress 


Finite  Elasticity: 

Consider  a  solid  body  with  an  initial  unde  formed  configuration  B  and  a 
deformed  configuration  b.  Let  the  coordinates  of  a  point  P  in  B  be  ^ ,  with 
primary  base  vectors  G^.  The  material  particle  at  P  in  B  is  assumed  to  have 
moved  to  the  point  p  in  *b'.  Let  *b*  be  defined  by  another  set  of  arbitrary 
curvilinear  spatial  coordinates  r|L,  or  by  convected  coordinates  £"*,  with  base 
vectors  and  ^  respectively.  Let  the  vector  of  displacement  of  the  particle 
from  P  to  p  be  u  =  u^G .  The  deformation  gradient  tensor  F  is: 

l  -  Ff)£1SJ  ■ 

-14%SJ  (i) 

I 

wherein,  use  has  been  made  of  the  dyadic  notation.  In  terms  of  jj,  we  can  write 
F  =  FLrg/  =  (5^+uVk)  G/  -  I(P)  +  e  (2) 

where  I(P)  is  the  identity  (metric)  tensor  at  P,  and  (  ) ,  denotes  a  covariant 

*«•  K. 

derivative  w.r.t.  t,  at  P.  We  see  that  the  displacement  gradient  tensor  e  is: 

5  -  J&f-  <« 

Let  r  be  the  tensor  of  true  or  Cauchy  stress  at  p,  which  measures  tractions  on 
an  oriented  surface  (dan)  at  p.  If  the  image  of  (dan)  in  the  undeformed 
configuration  is  (dAN)  at  P,  we  define  the  first  and  second  Piola-Kirchhof f 
stress  tensors,  denoted  here  by  t  and  s  respectively,  as: 

(dAN)  .  t  =  dan  .  t  (4) 

(dAN)  .  s  =  (dan)  .  (t.F_T)  (5) 

-1  T 

where  (  )  denotes  an  inverse  of  a  tensor,  and  (  )  denotes  a  transpose. 


Using  the  geometrical  relation 


(dan)  =  .1  (daN)  .  F 


(6) 


where  J  is  the  lbsolute  determinant  of  F  (or  the  determinant  of  the  matrix 

[FL‘])«  Eqs.  (4)  and  (5)  lead  to: 

#  K. 

t  =  JF"1  .  T  =  F_1  .  o  (7) 

-1  -T  -1  -T  -T 

s-JF.T.P-P.g.F-t.F  (8) 

where,  g  =  Jt  is  the  Kirchhoff  stress. 

Noting  that: 


-1XJ. 


(F  ).*GjS  - 


1C 

■Bn1 


Gri  =  G.j> 


(9) 


and 

i.  j  K*.  L 
~  =  T.j8i8  =  t.l*%8 

we  have : 


(10) 


J  — x*"G  & 

an1  -n-J& 


T  K*.,,  L 


(11a) 


and 


8  =  J  “  t™1'  ^  G  G  .  (lib) 

.  m  .  n  — J— L 

3n  3n 

Let  W  be  the  strain-energy  density  measured  per  unit  volume  at  P  in  B.  Then  it 
is  known  that : 


3W  _  3W  T 
3F  3e  ~  - 


(12) 


and 


3W  =  9  3W  =  s 
3^  ~  “  3C  ' 

T  1 

wliere  C  =  F  .F  and  y  =  —  [C—  I ( P )  ]  is  the  Green-Lagrange  strain 
equations  of  linear  (LMB)  and  angular  (AMB)  momentum  balance 
(LMB)  :  7p  .  t  +  oo(f-a)  *  0 


(13) 

tensor.  The 
conditions  are: 


or  Vp  .  [s.FT]  +oQ(f-a)  =  0 


(14a, b) 
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(AMB)  :  F  .  t  =  tT  .  FT 


or  s  =  a  . 


(15a, b) 


Compatibilicv : 


F  -  I(P)  +  e  •  (6^+u^)  -L-K 


lr  ,  '  M  ,  „L_K 

or  I  =  2[UL,K  +  “K.L  +  ^LUM,K]  ^  ’ 


Boundary  Conditions: 


N  .  t  =  N  .  (s . F  )  =  t  at  S 


u  =  u  at  S  . 
—  —  u 


(16a, b) 


(17a) 

(17b) 


For  our  present  purposes,  we  restate  Eq.  (17b)  in  a  bit  more  ambiguous,  but 


still  physically  meaningful  form  as: 


F  =  F  at  Sp. 


(17c) 


The  AMB  condition,  Eq.  (15a)  is  automatically  satisfied  provided  W  is  a 


frame-indifferent  function  of  F  [i.e.,  for  instance,  when  W  is  a  function 

of  F  only  through  UT(F)  =  W(C)  =  W(F^.F)]  and  t  is  defined  through  Eq.  (12). 

In  Eq .  (14),  p  is  mass/unit  volume  at  P  in  B,  f  are  arbitrary  body  forces 
o  — 

per  unit  mass,  and  a  is  the  absolute  acceleration  of  the  material  particle 

L  2  2  L  L  — 

(a  =  a  C  =  (d  u./dt  )  G  5  ii  G  ].  Likewise,  in  Eq .  (17a),  _t  are  prescribed 

L  Li  L 

tractions  per  unit  undeformed  area  S  of  B.  If  t,  s,  F,  e,  ii  are  fields  that 
satisfy  Eqs.  (12  through  17),  we  now  show  that  for  any  close  volume  V,  which 
is  free  from  anv  singularities  or  other  defects,  in  the  configuration  B,  the 
following  ident  it  ies*  are  valid: 

0  =  (t.F)-Po(_f-a).  F]dv  +  J  (N.t-t).Fds  +J  N.t.(F-F)ds  (ll 


0  =J  f  (s.C)-Po(f-a).  F|dV  +  f  (N.s.Fr-t).Fds  +  J  N.s.F 


i.F  . (F-F)dS 


*  These  represent  a  generalization  of  the  conservation  law  given  in  [3],  :o  the 

case  when  body  forces,  inertia,  and  arbitrary  crack-face  conditions  are  considered. 


Ill  the  above,  we  made  the  assumption  that  the  material  is  homogeneous,  i.e., 
W  is  a  function  of  the  location  P  in  B  only  by  virtue  of  the  fact  that  F  is, 
in  general,  a  function  of  P.  The  proof  of  Eq.  (18)  is  evident  from  the  fol¬ 
lowing  : 


V  =  CL  (-^~) 
P  ar 


r  t 

-  G  it  : 


•l]  ■  ?  [' 


=  GL  t 


:  (fkj>’l^J_ 

4 

>  4™<v 

•J 


where  (  ),^  denotes  covariant  differentiation  w.r.t.  £  .  And, 


Ip  •  <E'f>  -  ■  -rlfM\i&£K] 


ar 


=  ^  -  ,m  s.  rL 

(  FNK  '  (t  fn4 *m— 


MN 


MN„ 


But 


^  ’MFNL+t  FN1,M  — 


FNL,M  ~  ^NL^N.L^’M  ^N.L^’M  ^N.M^’L 


(GNM+UN’,M)  ’L  ^NM^’L 


(20) 


(21) 


(22) 


since,  the  metric  tensor  behaves  as  a  constant  during  covariant  differentiation. 
Thus,  equation  (21)  can  be  written  as: 

-p  •  (E-E>  ■  (t”‘-r,Fm.+tMNF»M.L)eL 


-  (Vp.t)  .  F  +  [tT:(F),L]GL.  (23) 

Upon  using  F.qs.  (20),  (23),  (14a),  (17a),  and  (17c),  the  validity  of  Eq .  (18) 

T  T 

is  immediately  evident.  When  it  is  noted  that  t  =  s.F  ,  and  F  .  F=C,  the  proof 
of  F.q.  (19)  is  apparent.  An  independent  proof,  follows  from: 

ipW  =  GL[s  :  (y),L]  =  GLU  :  j((KT),l  •  F  +  4  ■  (H.J) 
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=  CL{s  :  [ ( F)T ,  .F)]}  =  C  {  (s.FT)  :  (FT)  }.  (24) 

v  L  v  i-* 

Since  s  is  a  symmetric  tensor.  Further 

1?  •  [?.C]  =  Vp  .  [s  .  FT  .  F]  -=  [£p  .  (s .  F?)  ]  .  F  +  (s.FT  :  (F)  >L  CL 

=  [Vp  .  (s.FT)]  .  F+  (s.F?)  :  (FT),lGL.  (25) 

From  (24),  (25),  (14b),  (17a)  and  (17c),  the  validity  of  Eq.  (19)  can  be  noted. 
Using  the  Green-Gauss  theorem,  we  can  write  Eqs.  (18)  and  (19)  in  the 


0  =  f  [NW 

J3V 

+  f  N 

J sF  " 

•  L 


[NW  -  N  .  (t.F) ]  ds  -  / 

V  ' 


-  (  O  (f- 

J  V  ° 


a)  .  Fdv  +  /  (N. t-t)  .  Fds 

J  S~ 


t  .  (F-F)  ds 


NW  -  N  .  (s.C) 


-  f  p  (f-a).FdV  +  f  (N.t-t)  .  Fds 

Jv  °  J 


L  a,  i  ■ 


(F-F)  ds 


aero  IV  is  he  surface  of  V,  and  we  assume  that 


3V  =  S  +  S„  +  S.  . 
t  F  l 


Tlius,  Eq .  (26)  can  be  written  as: 


[NW  -  N  .  (t.F) ]  ds 


-  f  P  (f-a)  . 

Jv  °  "  “ 


+  J  [NW  -  l  .  F)  ds  +  J  (NW-N.t.F)  ds.  (29) 

S  t  'S  F 

Eventlnm.-.h  Eq .  (29)  is  applicable  to  a  general  )-d  i  mensiona  1  case,  an  illus¬ 
tration  of  the  2-dimensional  case  is  given  in  Fig.  I,  wherein  S.,  S_,  and  S 

i  F  t 

are  depicted,  and  N  is  an  unit  'outward1  normal  to  IV  as  shown  also  in  Fig.  1. 
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Suppose  we  consider  the  cartesian  coordinate  system:  along  the  crack 

surface,  x^  normal  to  the  crack  face,  and  x^  along  the  crack  front,  and  con¬ 
sider  the  component  along  the  x^  direction  of  the  vector  identity,  Eq.  (29). 

Let  the  problem  be  also  a  special  case  in  which  (i)  the  crack-faces  are  free 
of  tractions  and  any  imposed  displacement  conditions,  (ii)  the  body  forces  f^ 
are  zero,  and  (iii)  the  inertia  is  negligible,  i.e.,  the  problem  is  one  of 
elastostatics.  In  this  case,  since  =  0  along  the  face,  we  have  from  Eq.  (29), 


y-  yJkFki> ds :  X  (Ni“-  ViW  ds- 


It  is  this  sense  of  path-independence  of  the  integral  on  T  ,  and  the  as- 

oLU 

sociated  physical  interpretation  of  the  integral,  that  were  essentially 
presented  by  [2].  However,  in  the  general  case,  i.e.,  the  case  in  which:  any 
of  the  conditions  (i)  -  (iii)  above  are  not  satisfied,  and,  in  addition, 
the  compon  its  in  and  directions  of  the  vector  identity  (29)  are  also  of 
interest,  ve  need  to  r ’-examine  the  above  path-independence.  For  purposes 
of  understanding  the  ;.bove  general  case,  let  us  apply*the  result  in  Eq.  (29) 
to  a  vol  e  V-V  as  hown  in  Fig.  2a,  for  a  2-dimensional  problem  and  for  3-di¬ 
mensional  problem  in  Fig.  2b.  For  purposes  of  clarity  in  presentation,  let  us 
consider,  without  loss  of  generality,  the  2-dimensional  case.  In  view  of  Eq.  (29) 
we  then  have,  referring  to  Fig.  2a, 


[NW-N.(t.F)]  ds  + 


[NW-N. (t.F)l  ds 


-  f  p  (f-al 

Jv-v 


)  .  Fdv 


+  f  (NW-t.Fj  ds  +  /  [NW-N.t.F]  ds 

Jv  -'r” 

w  1  /.  c  1  /  r 


f  (NW-t.F]  ds  +  /  [NW-N.t.F]  ds  =  0 

Jr'  ~  Jr" 

»  1  I  -1  1  \  r* 


*  Note  that  the  divergence  theorem  cannot ,  in  general,  be  applied  when  the  volume 
integral  contains  a  non- integral’ le  singularity.  Thus,  when  volume  V  includes  the 
crack-tip,  referring  to  Eq.  (18),  q  W  mav  be  of  order  r”^  near  the  crack-tip  (front); 
and  hence  non- integrab 1 e .  Thus  the  divergence  theorem  cannot  be  applied  to  (Eq.  18)j 
in  the  case  of  V,  but  onlv  in  the  case  of  V-V  in  the  limit  ] 


where  r,,,  is  the  contour  with  a  unit  normal  directed  inwards  (into)  V-V 
obi  e 

as  in  Fig.  2a.  and  F^  are  portions  of  where  tractions  and  displace¬ 

ments,  respectively,  are  applied.  Similar  definitions  apply  to  r'^  and  F^* 
Since  the  crack  is  mathematically  a  surface  of  discontinuity,  it  is  seen  that 
N  ,  which  is  a  unit  outward  normal  to  is  equal  to  the  negative  of  N+, 
which  is  the  outward  normal  to  F^’  We  now  write: 

fv-V  P(A'?dV  =  fv-v  l(J^V  ~  ioM  dV 

E  E 

*  x_v  I!?  dv*  (32) 

G 

Using  the  divergence  theorem,  we  further  write: 

i 

J  Vp  (^P^u.cldv  =  -  J  N  (^p^-iOds  -  f  N  (|poa-u) 

V_Ve  °"~  ^r561  Jt21  4 

-J  I+[(|  Pau.u)+- (-|pou.u)_]ds.  (33) 

*^ri2 

Using  Eqs.  (32)  and  (33)  Li  Eq.  (31),  we  write: 


f  [N(W-T)  -  N  .  (t .  F)  1  ds  -  f  p  f  .  Fdv  +  f  jTT  (p  u.F)l 

J r„.  '  ~  Jv-v  '  Jv-v  ,ac  °  '  I 

234  g  g 

f  N+  [W+-W~)  -  (T+-T~)]  ds  -  f  [(t.F)+  +  (t.F)"]  ds 
*'ri2  r  12 

-  f  N+  .[ (t.F)+  -  (t.F)-]  ds 
Jr" 

1  12 

=  f  [N(W-T)  -  N  .  (t.F)]  ds.  (34) 

Jr 

c 

Note  that  in  the  above,  considering,  without  loss  of  generality,  the  segment 
12  to  be  of  the  same  length  as  the  segment  45,  the  integrals  on  and 

in  Eq.  (31)  have  been  combined  into  a  single  integral  on  alone,  with 


appropriate  sign  changes  in  the  unit  normal  being  accounted  for.  Further, 
the  notation  for  the  kinetic  energy  density,  T  =  (l/2)p^i.ii  has  been  employed: 
and  r  is  now  equal  to  with  a  unit  normal  acting  inwards  (into)  (V-V£) 

as  shown  in  Fig.  2a. 

We  now  consider  Eq.  (34)  in  the  limit  when  e-*-0,  i.e.  the  volume  ,  and 

the  contour  V  shrink  to  zero.  First  consider  the  right  hand  side  of  Eq.  (34). 

e  ; 

Let  ~  ,  in  the  2-dimensional  problem,  be  a  circle  of  radius  e.  (In  3-dimensional 
e 

problems,  we  let  r  be  a  circular  cylindrical  surface  of  radius  c,  with  the 
axis  of  the  cylinder  being  the  crack  front).  Then, 


Lt 

c-0 


-  Lt 

e-0 


f  [N(U-T) 

J  r 

e 

+n 


/ 


[N(W-T) 


-n- 


-  n  . 


-  N  . 


(t.F)  ]  ds 


(t.F)  ]  , _ .Ed0  . 

-  -  (r-E) 


(35) 


It  is  seen  that  the  right-hand  side  of  * q.  (35),  for  a  sharp  crack,  would  vanish 
identically  unless  W  (and  by  dimensional'  considerations,  t.F),  and  T  all  have 
singularities  of  the  type  r  a  (with  r,0  being  pola: -coordinates  centered  at  the 
crack-tip)  where  a>l.  However,  since  1  he  total  s' rain  and  kinetic  energies  in 
a  small  core  region  near  the  crack-tip  must  be  finite,  it  is  seen  that  “must 
be  equal  to  1.  Since  W  is  a  nonlinear  function  of  F,  F  can  be  expected  to  have 

_o 

a  singularity  of  type  r  with  ft<l;  likewise  t  may  have  a  singularity  of  type 

r  ^ ,  such  that  j+i=l.  It  is  worth  noting  that  available  solutions  in  linear 

elas todynnmic  crack  propagation  indicate  that  the  material  time  derivative  ii 

-1/2 

(or  absolute  velocity  of  a  material  particle)  may  have  a  singularity  of  r 
On  the  other  hand,  if  the  crack  is  stationary,  even  when  the  dynamic  effects 

■fl/a 

are  accounted  for,  u  varies  as  r  “[8].  Thus,  in  general,  in  a  dynamic  crack 


propagation  problem,  the  term  on  R.H.S.  of  Eq.  (34)  is  non-zero  for  a  sharp 


crack,  even  in  the  limit  as  e-K) .  We  denote  this  non-zero  limit  as  the  vector 
J.  Thus,  from  Eqs .  (34)  and  (35)  we  have: 


f  (N(W-T)  -  N  .  ( t .  F) )  ds  +  Lt  I-  f  p  f.Fdv  +  f  ■—  (p  u.F)  dv 

Jr234  ~  “  '  '  c->o  '  Jv-ve  °"  '  Jv-vedt  ~ 

■f  f  N+  [(W+-W~)  -  (T+-t“)]  ds  -  J  [ (t • F)+  +  (£.?)■]  ds 
JV12  JV12 

~  j  N+- [(t.f)+  ~  (t.F)“]  dsj 
Jr'l2 

X 


Lt  I  [N(W-T)  -  N  .  (t . F) ]  ds 
e-K)  J F 


(36) 


=  J. 

Note  that  the  limit  e-K)  in  the  term  within  brackets  {  }  in  the  extreme  left 


hand  side  of  Eq.  (36)  implies  that:  a)  in  the  volume  integrals,  a  volume, 

however  small,  near  the  crack  tip  must  be  deleted,  and  (b)  in  the  crack  face 


integrals,  an  area,  however  small,  neai  the  crack-front  must  be  deleted. 

Now,  in  view  of  Eq.  (39)  it  is  j=  een  that  th. ■  integral  on  the  extreme 
left-hand-side  of  Eq .  (36)  is  path-in  iependent  t  (Note  however,  the  volume 
V-V£  change  for  each  path).  It  is  worth  noting  that  Eq.  (36)  implies  that 
the  path  234,  as  in  Fig.  2a,  begins  at  the  point  2  on  the  lower  flank  of  the 
crack,  and  ends  at  point  4  on  the  upper  flank,  such  that  points  2  and  4  are 
equidistant  from  the  crack-tip.  However,  this  is  just  for  convenience,  points 
2  and  4  need  not  be  equidistant;  in  which  case  the  integrals  on  in  Eq.  (36) 
must  be  "split  up"  into  integrals  on  both  r  ^  and  as  in  Eq.  (31). 

Evidently,  if  one  uses  Eq .  (27)  and  repeats  the  steps  in  Eqs.  (31-36), 
we  obtain  an  equivalent  representation  for  a  path-independent  integral: 


[ N(W-T)  -  N  .  s 


(2>+[) 1  ds  +  Lt  (-  f  of. 

K)  '  Jv-V  °~ 

€ 


Fdv 


*  Note  the  fundamental  difference 

in  Eqs.  (30,  39)  (which  are  due 

(36)  on  the  other.  In  Eq .  (36) 

evaluated  on  the  contours  f  ,  _  , , 

12  34  a 


in  the  notions  of  pa 
fo  [2,4|)  on  the  one 
it  implies  that  the 
or  r  *  /.  •  c.  (see 


th-independence  as  embodied 
hand,  and  the  present  Eq . 
integral  on  the  extreme  L.H.S 
Fig.  2a)  has  the  same  value. 


1 2  ’  3  *  4  ’  5 


/  ^  (p  u.F)  dv  +  /  N+  (W  -W  )  -  (T+-T  )]  ds  -  f  [(t.F) 

Jv-v  dt  °“  '  J r  ~  Jr' 

e  12  12 

(t.F)l  ds  -  J  N+.[(s.FT.F)+  -  (s.FT.F)”]  ds  |  =  J.  (37) 


Trivial  as  the  difference  may  be  from  a  theoretical  view  point,  Eq.  (37)  is  in  fact 

t 

more  convenient  from  a  computational  view  point,  to  calculate  J  from  a  far- 
field  contour,  than  Eq .  (36),  since  most  generally  available  computer  programs 
(finite-element!)  use  s  and  y  as  primitive  variables  [13 . 

If  we  consider  the  special  case  when  (i)  deformations  are  infinitesimal, 
and  thus  the  distinction  between  various  stress  measures  vanish,  i.e.,  t=o,  (ii)  the 
material  is  linear  elastic  and  homogeneous,  (iii)  the  material  is  under  dynamic 
equilibrium,  and  (iv)  the  geometries  of  the  solid  and  crack,  are  conveniently 
described  in  a  cartesian  system,  then  Eq.  (30  is  reduced,  at  any  time  t,  to: 


JK  ‘  X  (VW-T)  -  Vmj'V-J.K51  ds 

r2  34 

+  ’  /v-v  dV  +  [_v  ii  lpoV5JK+uJ,K)|dV| 

e  e 

+  f  Nj  [(w+-w“)  -  (T+-T-)]  ds  -  f 

•r  ^  \  O  ~  ^  1 


l<tJ(6JK+UJ,K)>  +  <tJ(6JK+UJ,K)>  ]  ds 


fr 


(tJLFLK)  ]  ds- 


12 


(38) 


This  should  be  contrasted  with  the  expression  for  linear  elasto-static  case, 
when  no  conditions  are  prescribed  on  the  crack-face,  given  in  [2,4]  denoted 
here  for  comparison  purposes,  as  J*: 


jr 


-i. 


(wVVmjuj,k]  ds- 


(39) 


34 


Thus,  in  the  present  J,  Eq.  (38),  discontinuities  of  W  and  T  across  the 
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crack- faces* a re  allowed  for,  in  addition  to  its  being  applicable  to  the  general 


case  as  discussed. 

It  is  interesting  to  note  however,  that  even  in  the  general  case  repre¬ 
sented  by  Eq.  (36),  the  following  identity  holds: 


f  Nr  .5,,,  +  Lt  (-  f  p  f  5  w 

Jr,,.  M  w  JK  e->o  l  Jv-v  0  J  JK 

234  e 

r  v  &  (poVjk)  dv  -  f,  <  (?j6jk)+  +  <Vjk>’  >ds) 

J  e  JV12 

fr 


=  Lt 

£-*■0  Jr 


!U6JKdS 


(40) 


since,  Eq.  (40)  is  nothing  other  than  the  global  linear  momentum  balance  con¬ 
dition  for  the  domain  V-V^  in  the  limit  e->-0.  Thus  we  may  reduce  Eq.  (36), 
through  a  slight  modification  to  the  definition  of  J,  as: 


f 


[N(W-T)  -  N  .  (t.e)]  ds  +  Lt 


234 


e-K) 


/ 

Jv-v 


p  f.edv 


f  j-  Cp  u.e)  dv  +  f  N+[(W+-W")  -  (T+-T  )]  ds 

Jv-V  ac  “  Jr., 

£  12 


f  +  (£•?)  1  ds  - 

f  N+.  [(t.e)+  -(t.e)-]  dsl 

^12 

Jr"  ~  ~  I 

7  12 

=  Lt 
e-K) 


t. 


[N(W-T)  -  N  .  (t.e)]  ds 


£  J. 


(41) 


Thus  Eq.  (41)  represents  a  slight  modification  to  the  definition  of  J_  in  as  much 
as  F  in  Eq.  (36)  is  now  replaced  by  e[e=F-T]. 

We  now  examine  the  physical  interpretation  of  J  as  presently  defined 


through  Eq.  (41).  For  this  purpose  let  us  consider  the  volume  V  at  time  t. 

c 

It  is  to  be  understood  that  V  is  a  small  region  at  the  crack-tip  (front)  with 

*  Note  that,  based  on  physical  cons iderat ions ,  we  may  assume  that  W+  and  W  along 
the  crack- face,  in  general,  may  have  integrable  singularities  near  the  crack- 
tip  (front).  An  example  wherein  (W^-W)  ^  0  at  the  crack-face  has  recently  been 
brought  to  the  author's  attention  [13]. 


tue  surface  r  .  Let  p  be  the  potential  of  external  forces  acting  on  V  ,  and 

£  £  C 

let  W*  and  T*  be  the  strain  energy  and  kinetic  energy,  respectively  of  .  It 
is  seen  that 


■  -  /v  !ivv  -  /s  v*i dB-/r  tr0  <42) 

e  te  e 

r 

where  S  is  the  crack  sutface(s)  enclosed  by  T  .  The  last  term  on  the  r.h.s. 

tlC  £ 

of  Eq.  (42)  is  the  work  of  tractions  exerted  on  by  the  surrounding  solid, 
and  these  tractions  are  dependent  on  the  displacement  field.  Likewise,  we  have 


(43) 


and 


T*  =*  f  Tdv  *  f  |pu  u 

e  Jv  Jv  2  1  1 


dV 


(44) 


where  u^  is  the  absolute  velocity  of  a  material  particle.  In  ?•  s.  (42-44)  car¬ 
tesian  coordinates  x^  have  been  used  for  simplicity.  Let  c^  be  the  cartesian 
coordinates  of  the  crack-tip .  Let  i;^(=x^-c^)  be  coordinates  ct  ntered  at  tie 
crack-tip  when  the  crack  is  of  "length"  c^.  It  is  seen  that  i  >.  general,  in  the 
immediate  vicinity  of  the  crack-tip,  i.e.,  in  V  ,  which  may  be  considered  to 
be  a  "process-zone",  we  may  assume: 


U1  =  Ui^r'k’Ck’^k’ =  tij  ^r’k,Ck’^k,t:^ 

W  =  W(^k » c^,  ) ; T  ~  T (4j, » ci, , i ). 


k’  k’  k’ 


(45) 


Thus  the  variables  u^,  t,  W,  and  T  may  depend  explicitly  on  the  crack  "length" 
c^,  crack-velocity  c^,  as  well  as  time  t.  It  is  known  that  in  the  vicinity  of 
a  sharp  crack-tip  (front),  the  displacements  u^  are  nonsingular  (even  though 
the  material  velocity  iV  may  be  singular),  while  W  and  T  are  singular:  and 
based  on  physical  considerations,  the  singularities  in  W  and  T  can  be  of  order 
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r  1  where  r,  9  are  polar  coordinates  centered  at  the  crack-tip.  Even  in  finite 
bodies  the  representations  (45)  may  be  considered  as  being  asymptotically  cor¬ 
rect  within  the  "process-zone"  irrespective  of  the  geometry  of  the  body  and 
loading.  We  may  write  at  time  t,  for  instance. 


( 


W* 


/  WUK’CK*£k’t)  dV 


(46) 


£ 

and  a  similar  expression  for  T*  of  v  . 

c 

Now,  let  us  assume  that  in  time  (dt)  the  crack  "advances"  by  dc^.  Let 

£k  be  once  again,  the  cartesian  coordinates  centered  at  the  new  crack-tip  such 

that  <;k=r  -dc  as  shown  in  Fig.  3.  Also,  without  loss  of  generality,  let  £=6dc^, 

8>1,  such  that  the  crack-tip  which  is  within  (say  a  circle  of  radius  e  centered 

at  the  crack-tip)  at  time  t,  will  not  lie  outside  the  fixec.  at  time  (t+dt) . 

Thus  we  are  considering  the  case  when  the  crack-tip  advances  by  "dck"  in  time 

'dt'  into  a  fixed  volume  V  .  Analogous  to  Eq.  (46),  we  may  then  write: 

-  e 


W* (t+dt)  = 


W(c^,ck+dck,ck+dck,t+dt)  dv. 


(47) 


£ 

We  note  that,  in  the  asymptotic  sense,  the  dependence  of  W  on  c, £  at  t+dt  for 
a  sharp  crack  is,  in  general,  of  the  same  functional  form  as  the  dependence 
of  W  on  'k  at  t.  Note  also  that  in  performing  the  integration  as  in  Eq.  (47), 
the  limits  of  integration  for  a  f ixcd  volume  V  in  terms  of  would,  of  course 
be  different  from  those  in  terms  of  'k  in  Eq .  (46)  (see  Figs.  3a  and  3b).  How¬ 
ever,  by  noting  that  f -dc,  ,  and  dr'=dr,.  ,  we  have,  by  what  amounts  to  a 

k  |c  k  k  k 

change  of  variables. 


W* ( t+dt )  =  /w[(,k-dck),(ck+dck).(ck+dck).(t+dt)l  dv. 


(48) 


Now  the  limits  of  integration  in  Eq.  (48),  would  be  the  same  as  those  in  Eq. 


(47).  Now,  we  may  expand  out  the  integrand  on  r.h.s.  of  Eq.  (48)  using  Taylor 


series  and  express  it  in  terms  of  W(^^,c^,c^, t) .  However,  it  is  recalled  that 
W  is  singular  (viz.,  r  ^  for  a  sharp  crack)  w.r.t.  as  4^"^.  where  as  the 
explicit  dependence  of  W  on  c^,  c^  and  t  is  in  general,  non-singular.  Thus, 
while  DW/3c^,  3W/3c^,  and  3W/3t  are  all  integrable,  the  partial  derivative 

t 

3W/3c^  would,  however,  be  non-integrable .  Thus  the  above  idea  of  changing 
variables  is  non-workable .  Thus,  we  rewrite  Eq.  (47)  as: 


W*(t+dt)  = 

c 


Jr  tut  .  .  ^  3W  .  ,  3W  ,.  ,  3W  . 

[W(c,'  ,c,_ ,c,  , t)  +  dc,_  +  ^-r-  dc,_  +  —  dt]  dv. 

V 


*k’  k’  k’ 


3  c, 


3c, 


3t 


(48) 


Now  consider  the  term: 


=  X  W(^’Ck,dk’t)  dv  ’  X  W(WVt)dv* 


(49) 


Since  the  functional  dependence  of  W  on  sk  and  £k  are  of  the  same  form,  we 
can  "subtract  out"  the  singularities  in  evaluating  I  of  (49)  as  sh jwn  in  Fig. 
3c.  From  Fig.  3c  it  is  apparent  that  the  term  I  of  (49)  is  given  by: 


I  = 


J  (WNkdcR)  ds. 


(50) 


e 

Note  the  negative  sign  on  the  r.h.s.  of  Eq.  (50)  is  due  to  the  definition  of 
the  "outward"  normal  to  the  contour  T  in  the  sense  indicated  in  Fig.  3c.  [i..e, 

a  contour  beginning  at  the  "bottom"  surface  of  crack  and  ending  on  the  "top" 
surface].  Let  us  now  define  the  derivative  [D(  ) / Dt ]  as  the  total  rate  of 
change  of  (  )  in  a  time  "dt"  due  to  crack  growth  by  "dck".  Thus,  for  instance, 


« 


DW* 

{-fr^lCdt  =  W*(t  +  dt)  -  W*(t). 

UC  C  G 


From  Eqs.  (48),  (49)  and  (50)  it  is  seen  that: 


DW* 


WN,  c  ds  + 
k  k 


(51) 


(52) 
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In  the  case  of  a  dynamically  propagating  crack,  the  kinetic  energy  T  may  also 
have  (as  is  known  [8],  for  instance,  in  linear  elastodynamics)  a  singularity. 
Thus,  following  essentially  the  same  arguments  as  in  Eqs.  (46-50),  we  have: 


DT*  r 

•  -  X 


™  J  ,  f  /  3T  .  ,  3T  ..  ,  3T.  , 

TN,  c  ds  +  I  (t —  c  +  rr-  c  +  — )  dv. 

k  k  jv  3ck  k  3ck  k  3t 


Finally,  at  times  t  and  (t+dt),  respectively,  the  displacements  are  given 

Du. 

L  c 

by  u^  and  [u^+(-jj^~)  dt].  From  the  F.quations: 

ui[tl  =  V'WV0 


u^ [ t+dt ]  =  ui(^k,ck+dck,ck+dck,t+dt) 


=  ui(4k-dck,ck+dck,xk+dck,t+dt) 


it  is  seen  that: 


Du. 

(— i)c  = 

V  Dt  } 


3u^  3u^  3u^  3u^ 

3^  dk  +  3^  Ck  +  36^  Ck  +  IT* 


Using  Eq.  (55),  it  is  then  seen  from  (42)  that: 

DO  /-Du  r  Du  /•Du. 

<dt>c  -  -/v  !!  <-5r'° -  X  J1  ("5r,c““'X  h  (-ST>Cds  '  6) 

c  tec 

where  t.  are  tractions  corresponding  to  u.  at  T  at  t.  We  will  now  consider 
l  i  e 

the  "energy  release  rate"  to  the  crack-tip,  as  measured  in  the  process-zone, 

c 

.  This  energy  release  rate,  denoted  here  by  (DE^/Dt)  ,  is  given,  from  the 
energy  balance  within  V^,  by: 

DE  Dt|i  DW  DT 

(sr)C  -  -  (-dt)C  -  [(-dt)C  +  (3T)]-  (57) 

Upon  using  (52),  (53)  and  (56),  we  write  (57)  as: 

DE  |  f  Ju.  r  3u.  r  3u 

<-  V  .  £k  I  /  l(w+T)Nt -.aid.-  /  f.  ja  dv  _  /  _ 
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■  Cfc  l/v 

" tk  I L 


3W  _3T_ 

V  3ck  3ck 
c 


( aw  _<rr_ 
v  3‘k  3dk 

e 


y*  3u 

— i)  dv  - 
1  3ck 

JT  li  3cj 

£ 

8u. 

3u 


£ 

ir  .3»ji  ?  aV  r  3ui  r  3ui  j 

-  \L  («  +  «  -  ?i  ir> dv  -  X  i  it  ds  ■  X  ^  I 


te 


te 


3u , 


'v .  Jr 

e  e 


te 


e  (I)  +  (II)  +  (III)  +  (IV)  (58) 

where  (I),  (II),  (III)  and  (IV)  identify,  respectively  the  first  through  fourth 

terms  each  enclosed  within  {  },  on  the  r.h.s.  of  Eq.  (58).  Suppose  now  that 

for  a  dynamically  propagating  crack,  W,  T,  and  t„  au^/ax^  [note  that  3(  )/3x^ 

E  3(  )/3^i]  all  have  singularities  of  the  type  r  On  the  other  hand,  the 

explicit  dependence  of  these  quantities  on  c^,  c^,  and  t  are,  in  general,  such 

that  their  partial  derivatives  w.r.t.  c^,  c^,  and  t  are  also  singular,  but  the 

-1 


singularity  is  still  of  order  r 


Suppose  we  consider  a  2-dimensional  problem 


wherein  V  is  a  circular  domain  of  radius  e(-8dc  )  centered  at  the  crack  tip: 
e  k 


Ir 

€ 

f 


+n 

)  ds  —  I  ( 


)dv 


/ 

-n 

£  ' 

■  If 


)ed0 


n 

e  +F 


) rdrdO . 


(59a) 


(59b) 


o  -n 


Thus,  for  a  2-dimensional  problem  (the  same  argument  carries  over  to  3-dimen¬ 
sional  case)  for  instance,  the  term  (I)  of  Eq.  (58)  becomes: 

+n  „  c  +n 


(i)  = 


I  **  3u .  /•  3u , 

‘  M  /  «««>\  -  h  '  J  J  h  3T  'dtd8 

-n  *•  ’  o  -n 


t. 

-/ 


-  ui  -  3ui 

(ti  TT*  +  (ti  I  dr 

1  J  k  -[,  1  -  'k  +11 


(60) 


S  is  defined  by  a  =  +fl .  Supposing  that  f.  and  t  are  non-singular 
te  -  j.  ^ 


where 


1 


near  the  crack-tip,  while  W,  T  and  t  u  ,  behave  as  r  ,  it  is  seen  that  the 

X  1 ,  K 

first  term  on  the  r.h.s.  of  (60)  is  independent  of  c,  while  the  other  terms  are 
0(e).  Likewise,  it  is  observed  that  terms  (II),  (III)  and  (IV)  of  Eq.  (58a) 
are  all  of  0(e).  Thus,  in  general,  we  have  from  (58), 


^c-M/r 


[ (W+T)N  ~  t,  — ^  ]  ds  +  0(e). 
;  k  i  3^  J 


Hence,  in  the  limit,  we  may  take  for  an  arbitrarily  small  volume  near  the  crack- 


DE  r  3u.  . 

<-Df)C  -  Jk  1  J„  !<U+T)Nk  -  dS|- 


r  3u . 

/  [<W>Nk  -  tt  ds. 


(62a) 


Comparing  Eqs .  (62a)  and  (41)  it  is  seen  that: 


<-5T>C  ‘  Jk  Jk  +  2 


f  TVs|. 


Thus,  in  the  case  of  propagating  cracks  in  elastodynamic  fields,  Che  path  in¬ 
dependent  integral  of  Eqs.  (41),  and  hence  Eq .  (38),  does  not  have  the 
physical  meaning  of  energy  release  rate  per  unit  crack-growth.  However, 
if  we  define  the  "Lagrangean"  of  the  domain  V  as: 

L  =  -  il>  -  W*  +  T*  (64) 

e  e  c  c 

it  can  be  immediately  seen  from  Eq .  (57),  and  the  development  that  follows 
thereon,  that: 

DL  DE  DT* 

<dT)C  ■  (T5T'C  +  2<sr>c 


ik  /  I(u-T)Hk  -  t.  5^1 

k 


(h4b  ; 


•)  > 


Thus  the  path-independent  integral  of  (41),  and  hence  (38),  while  not  represen¬ 
ting  an  "energy  release",  does  still  have  a  physically  meaningful  interpretation 
of  rate  change  of  Lagrangean  per  unit  crack  growth. 

On  the  other  hand,  for  stationary  cracks  in  dynamic  elastic  fields  T  is 
nonsingular,  while,  in  elasto-statics  T  is  negligible  altogether.  In  these  cases, 

it  can  be  seen  that: 

DE 

<3T)C  dt  "  dVk  (65) 


and  thus,  does  have  the  physical  meaning  that  it  is  the  energy-release  rate 
per  unit  movement  of  the  crack-tip  in  x^  direction. 

In  the  linear-elastic  case,  an  expression  for  energy-release  rate  for  a 
crack  extending  in  an  e lasto-dynamic  field  was  derived  by  Freund  [7]  to  be: 


F  =  Lim  j  (o . .n.U.  +  4  a. .U.  .v  +  4  pU.U.v  )ds. 

^0  Ja  J  1  2  ij  i,j  n  2  i  i  n 


(66) 


s-*-0  Js 

where,  in  the  notation  of  [7J:  s  is  a  ''small’’  loop  that  moves  along  with  the 

crack- tip;  o  is  the  stress  tensor;  lh  are  displacements;  are  velocities, 

n^  are  direction  cosines  of  a  unit  normal  to  s  pointing  away  from  the  crack- 

tip,  and  v  is  the  "component  of  velocity  of  a  point  on  s  in  the  direction  of 
n 

n^  (if  the  crack-tip  is  moving  in  the  x^  -  direction  with  instantaneous  speed 
v,  then  vn=vnp-  1°  deriving  the  above  expression  (66),  a  small  loop  of 
unchanging  geometry  was  supposed  to  be  moving  with  the  crack-tip  and  a  global 
energy  rate  balance  was  employed  [7], 


It  is  seen  that  the  presently  derived  expression  for  energy  release  rate 
viz.,  Fq.  (62a)  differs  from  that  derived  in  [7],  viz.,  (66).  As  demonstrated 
earlier,  the  present  Eq.  (62a)  can  be  readilv  reduced  to  the  well-known  results 
for  energy  release  rates  for  elasto-stat ic  crack  problem  as  well  as  for  the 


elasto-dynamic  problem  when  the  crack  is  stationary.  On  the  other  hand  Eq. 
(66)  of  (7)  is  not  readilv  reducible  to  the  above  cases,  without  invoking 
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§ 

i 


certain  other  assumptions. 

We  now  consider  finite-strain  inelasticity,  but  restrict  our  attention 
to  the  two  cases:  (i)  finite  strain,  classical  rate-independent  elastoplas- 
ticity,  and  (ii)  finite  strain  rate-dependent  elasto-viscoplasticity 
including  creep  as  a  special  case. 

Finite-Strain,  Rate  Independent,  Classical  Elasto-Plasticity : 

Because  of  the  author's  own  interest  in  computational  mechanics,  the 
following  is  presented  in  a  fashion  that  is  directly  amenable  to  computations 
based  on,  say,  the  finite  element  method.  However,  the  basic  development 
itself  is  divorced  from  any  "computational"  overtones,  lest  they  may  be  deemed 
"bad". 

At  time  t,  the  material  particle  P  of  B  is  located  at p  in  b.  At  time 
t  +  At,  let  the  same  material  particle  move  to  p^.  As  discussed  before,  the 
spatial  coo»-  'inates  in  b  are  nm,  and  the  convected  coordinates  are  ^ .  Let  the 


vector  pp^  bo  Au. 

Au  =  .vu  C  =  ,Au  g  =  Au  g  .  (67) 

—  — J  -“m 

We  define  ne  incremei  tal-displacement-gradient  tensor  Ae  as: 

,  .  m  n 

Ae  =  Au  g  g 
min-2 

=  (Ae  +Aw  )  gV  (68) 

mn  inn 

where  Ae  and  Aw  are,  respectively,  the  incremental  strain  and  incremental  spin 
referred  to  the  current  configuration  at  t  and, 


L ,  ,  n  in  I  .  .  m  n 

U:  =  r(Au  +Au  )  £  g,  ;  Aw  =  (Au  ~Au  )  &  £ 

2  m;n  n ;  m  2  2  -  2  m;n  n ;  m  2  2 


In  the  above  and  in  the  preceding,  (  );  and  (  );  refer  to  covariant  deriva- 

n  .1 

tives(w.r.L.  r|n  and  ^  respectively  it  p,  using  the  metric  g  ,  and  g,„  res- 

mn  Jit 

pectivciy) in  the  current  configuration.  Let  the  Cauchy  stress  and  Kirchhoff 
stresses,  t  and  j  respectively,  at  time  t  be  represented  as: 


mn  J*K* 

T-  =  T  Vn  =  T  %% 

T  T  mn  T  J*K*  , 

=  Jr~  =  Jx~  =  jT  W  (70) 

Let  Ao  be  the  total,  or  substantial,  or  material  increment  of  Kirchhoff  stress, 
such  that. 


Ao  =  gAt  = 


(71) 


where  a  is  the  material  derivative  (for  a  fixed  material  particle,  ^  = 
const,),  of  g.  For  an  objective  stress-rate,  to  be  used  in  the  incremental 
constitutive  relation  for  an  elastic-plastic  solid,  we  take,  following  Hill 
[L4] ,  the  co-rotational  increment  of  o  (which  is  often  referred  to  as  the 
Zaremba,  or  the  Jaumann,  or  the  rigid-body  rate),  which  is  denoted  here  by 

A 

Ao,  It  is  well-known  (see  for  instance  fl_4] )  that 

Aa  =  Aa  +  a  .  Aw  -  Aw  .  a.  (72) 

A  constitutive  law  for  the  rate-theory  of  plasticity,  as  suggested  by  Hill 
[14]  is: 


A 

Ao 


3AV 

3Ae 


(73) 


where 


AV 


1  mn  pq 
—  L  Ac  Ac 

2  mnpq 


kL  2 

(i)(Xk£AE:  }  ' 


(74) 


Eq.  (74)  leads  to  a  bilinear  relation.  In  Eq.  (74):  L  is  a  tensor  of  in- 

mnpq 

stantaneous  elastic  modulii,  which  is  +  ve  definite  and  symmetric  under  mn->~t-pq 
interchange;  et=l,  or  0  according  to  whether  is  +  ve  or  -  ve,  X 

mn 

is  a  tensor  normal  to  the  interface  between  elastic  and  plastic  domains  in  the 
strain  space,  and  g  is  a  scalar  related  to  a  measure  of  rate  of  hardening  due 
to  plastic  deformation.  For  classical  isotropically  hardening  materials, 
the  above  constitutive  Law,  which  has  been  used  by  several  authors  in  the 
past  few  years,  as  discussed  inf  15],  becomes: 
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(75) 


12cxy^ (Ac :a ' )a ' 

Aa  =  2yAe  +  X(Ae:I)I  -  - 

(a':a')[6y+2  0F/3Wp)] 

where  X  and  p  are  Lame's  constants;  o'  =  a  -  -|(a:I)I  is  the  deviatoric 

% 

KIrchhoff  stress  at  p,  and  the  yield-surface  is  represented  by:  F  =  [3J0(a’)] 

Z  - 

-  F  (Wp)  =  0;  WP  =  /a:AePdt  and  J  (o’)  =  (1/2) (o ' :a ' ) . 
o  ,  1 

For  purposes  of  the  ensuing  discussion,  it  should  be  noted  that  the  con¬ 
stant  u  (which  is  equal  to  1  or  0)  in  Eqs.  (74,  75)  is  a  function  of  the  spatial 
coordinates  n"1  in  b.  Thus,  a  generic  point  p  in  b  may  be  experiencing  loading 
(a=l)  or  unloading  (a=0) .  Thus  at  time  t,  a  depends  on  the  location  of  p  in  b. 

Let  (x+At)  represent  the  first  Piola-Kirchhof f  stress  at  p^  at  time  (t+At), 
as  referred  to  the  configuration  b  at  time  t.  Then  it  is  seen  [15]  that: 

At  =  (Ao-Ae.o)/J  (76) 

=  (Aa-Ae .a-a . Aw) /J  (77) 

where  J=(p  /p  )  the  ra  io  of  mass  densities  in  B  and  b  respectively.  Use  has 
o  t 

been  made  of  Eq.  (72)  in  obtaining  (77)  from  (76).  Hence,  in  view  of  Eq.  (73) 
it  is  seen  [151  that'  a  potential  Aj  exists  such  that: 


and 


3AU 

T 

3A  e 

AU  = 


=  At 


AV 


1  T 

(Ac. Ac)  +  ~  q'  (Ae  .Ae) 


(78) 

(79) 


Note  AU  =  AU[a(p)],  i.e.,  AU  is  a  function  of  tx(=0  or  1)  which  is  in  turn 
a  function  of  p. 

Let  a  be  the  absolute  acceleration  (i.e.,  the  material  derivative  of  the 
velocity  vector),  of  the  material  particle  at  (t+At).  The  acceleration  a  can 
be  expressed  in  terms  of  spatial  bases  at  t,  or  sometimes  more  conveniently 
(from  a  computational  viewpoint)  in  the  fixed  bases  U  it  P  at  t=0.  (see,  for 
instance,  [15].)  The  field  equations  at  t+At  can  now  be  written  ns: 
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(80) 


(LHB)  V  .  [t+At]  +  P  (f-a)  =  0. 
— t  '  t  ~ 


(AMB) 

Ae  . 

T  T 

x  +  At  =  At  4-  t  .  Ae 

(81) 

II 

< 

3AU/3Ae^;  Ae  =  Au  gmgn 
m'.n^  ■“ 

(82,83) 

[x+At]  =  t  at  S 

(84) 

and 

Ae  = 

Ae  at  S  .  - 

(85) 

Note  that  V  =  gm(3/3nm)  is  the  gradient  operator  at  p  at  time  t.  In  the 
— t 

above,  t^  are  prescribed  tractions  at  (t+At)  as  measured  per  unit  area  of  surface 
S  at  time  t.  It  should  be  noted  that  when  the  potential  AU  for  At  is  as  defined 
in  Eq.  (79),  the  AMB  Eq.  (81)  is  automatically  embedded  in  the  structure  of  AU, 
as  discussed  in  [15]. 

« 

Now  we  show  that  for  a  closed  volume  v  (at  time  t)  wherein  a=constant 
(i.e.,  a  volume  which  is  everywhere  experiencing  loading  or  everywhere  exper¬ 
iencing  unloading),  and  which  is  free  fron.  any  singularities  or  other  defects, 
the  following  vector  integral  is  zero. 


X  I- 


(aiAe+AU)  -  (V  x)  :  Ae  -  V  .  [(x+A1.)  .  Ae] 


-  p  (_f-a)  .  Ae}  dv  + 


~|  dV  +  f, 


(x+At)  -  _t]  .  Aeds 


+x  - 


(x+At)  .  (Ae-Ae)  ds. 


(86) 


Note  again  that,  in  general,  AU  =  AUla(p)]  but,  within  the  considered  v^, 
a  =  constant  =  1  or  =  0.  The  proof  of  the  above  statement  is  based  on  the 
results : 

V  [ t : Au 1  =  (V  t)  :  Ae  +  [t  :  (  o)  ]  gm  (87) 

-t  -  — t-  '  ~  ;m“ 


v  AU  =  [AtT  :  (Ae)  ]  Rn. 
-  -  ;  m 


(88) 
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Since , 


V^a(p)  =  0  for  considered  (89) 

V  .  [(x+At)  .  Ae]  =  [V  .  (x+At)]  .  Ae  +  [(x+At)T  :  (Ae);  ]  g™.  (90) 

From  Eqs .  (87-89),  (80),  and  (82-85),  the  validity  of  Eq.  (86)  is  immediately 
evident.  , 

For  instance,  in  a  2-dimensional  dynamic  elasto-plastic  problem,  con¬ 


sider  two  path£*  T*  and  F*  surrounding  the  crack-tip  as  shown  in  Fig.  1  . 

We  note  that  within  the  region  (v^)^  enclosed  by  F*,  there  can  be  regions 

of  both  plastic  loading  and  elastic  unloading,  i.e.,  a  can  be  either  1  or 

0  at  various  points  within  (v^)f.  However,  within  the  region  (v^-v^)^, 

the  region  between  the  paths  T*  and  F*,  there  can  be  either  only  loading  or 

only  unloading,  i.e.,  within  the  region  (v^v.^,  either  <*=1 only  or  a=0  only.  In  this 

case,  we  have  in  view  of  Eq.  (86): 

0  =  f  ■  V(x:Ae+AU)  -  (V  t_)  :  Ae  -  V  .  [(x+At)  .  Ae] 

*'(v2-v1).  1 

-  p  ( f-a)  •  Ae  !  !v+  {  [  t i  .  (  r+At)-_t] .  Ae  ds+  f  n  .  (x+At) .  (Ae-Ae)ds  (91) 

t-  ')  l  2  ’  "*■  '  /  ' 

■'Cc 


<v; 


(se>l 


Recall  that  in  the  case  of  an  elastic  material,  the  basic  result  of  vanishing 
of  a  certain  integral  in  a  volume  free  from  singularities,  i.e.,  Eq.  (26),  was 
applied  to  a  region  V-V  ,  and  the  integral  vector  J  was  defined  as  the  non¬ 
vanishing  limit  of  a  contour  integral  on  F  in  the  limit  c-*0,  as  in  Eqs.  (36) 
or  (Al).  Note  that  Eq .  (91)  of  the  elasto-plastic  case  is  entirely  analogous 
to  Eq.  (26)  of  the  elastic-case.  However,  in  order  for  Eq .  (91)  to  apply  to 

the  region  V-V  in  the  limit  c-»0,  the  entire  region  (V-V  limit  e-O)  must  under- 

e  c 

go  either  plastic  loading  only,  or  elastic  unloading  (or  remain  purely  elastic) 
only.  Such  a  situation  is,  in  general,  unlikely  in  the  case  of  (growing)  cracks 


in  ductile  materials,  for  a  general  domain  V-V  .  If  indeed  these  conditions 

_  f. 

*For  convenience,  ant|  i  BCD  of  pig.  I  are  now  renamed  r^andr$,  respectively 


are  met,  and  hence  Eq .  (91)  Is  applicable,  in  the  domain  V-V^  (Lt  c-* 0)  ,  then 
results  entirely  analogous  to  Eq.  (36)  or  (41)  can  be  derived,  and  a  vector 
which  is  the  limiting  value  of  an  integral  on  in  the  limit  £-* 0  can  be  iden 
tified.  However,  in  general,  let  us  consider  the  case  when  the  conditions 
for  applicability  of  Eq.  (91)  are  not  met  in  V-V^  (Lt  e-*-0)  i.e.,  7^_a(p)?*0  in 
V-V  .  Referring  to  Fig.  2,  it  is  seen,  in  the  limit  e-*0. 


Lt  I  J  |v  (t:Ae+AU)  -  (V  t):Ae  -  V  .  [(x+At)  .  Ae] 
e-0  [Jv-V£  *~t  -t-  -  -t 

-  p  (f-a)  .  Ael  dv  +  f  ,  [n  •  (r+At)  -  £]  . 

c  ' 1  J(s  y 

t  e 

f  2  H 


Aeds 


t  0  =  R  (say) 


+  J  „  n  .  (x+At)  .  (Ae-Ae)  ds 
'e'  c 

where  R  can  be  seen  to  be  an  integr  1  over  V-V£  of  terms  involving  ^.a(p) 
Using  the  divergence  theorem  we  can  : ewrite  Eq.  (92),  for  instance,  for  a 
2-dimensional  domain  indicated  in  Fi.r-.2a, 


(92) 


f  [n  (x.-Ae+AU)  -  n  .  <  (fit)  .  Ae>]  ds  +  Lt  if  [-V  x:Ae 
J r„.w  ~  '  '  e-+0  l  JV-V£ 

:  ]  dv  +  f  n+  [(x:Ae+AU)+  -  (x:Ae+AU)  )  ds 
’  ‘  12 

J[(t.Ae)+  +  (t_.Ae)”]  ds  -  f  n+  [<(i+At)  .  Ae>+  -  <(r+At).  Ae>“]ds 
ri2  ri2 

=  R  +  Lt  r  [n  (x:Ae+AU)  -  n.  .  <  (x+At)  .  Aep]  ds. 

t:-K)  Jr 


234 


Pt(£-a)  .  Ae; 


AT  (say).  (93) 

It  is  noted  that  comments,  concerning  ‘  ^12'  ^12  eCc'*  essentiaHy  similar 

to  those  made  in  connection  with  Eq.  (36)  apply  in  the  case  of  (93)  as  well. 

In  view  of  Eq.  (91)  it  is  seen  that  the  integral  on  the  extreme  L.H.S. 
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of  Eq.  (93)is  "path-independent"  (however  the  enclosed  volumes  V-V  will  be  dif- 

c 

ferent)  for  any  two  paths  and  r2’3,4’  between  w^ich  there  is  entirely  loading 

or  entirely  unloading  i.e.,  a=l  or  ct*0  everywhere  between  the  two  paths.  Speci¬ 

fically,  consider  the  situation  when:  (i)T  ^  encloses  a  volume  (including 
the  crack-tip)  where  there  may  be  both  loading  and  unloading  taking  place  and 
(ii)  ^2' 3 '4'  encl°ses  a  volume  v2  such  that  in  v2~v^  also»  hoth  loading  and  un¬ 
loading  are  occuring.  Then  the  vector  integrals  (AT)^  (involving  ^234  on  t*ie 

extreme  L.H.S.  of  Eq.  (93))  and  (AT)2  (involving  r2' 3'4'^  not  be  equal. 

On  the  other  hand,  in  the  situation  when  (1)^234  encloses  v^  wherein  there 
may  be  both  loading  and  unloading  occuring  and  (ii)!^,^,^,  encloses  v^  such  that 
in  (v2_v  ) ,  either  only  loading  or  only  unloading  (or  remain  elastic)  is  taking 
place,  the  integrals  AT^  and  AT2  are  equal: 

It  is  noted  that  A1  integral  as  evaluated  on  ^234  fromEq.(93)  "measures "not  only 

the  severity  of  condition^  near  the  crack-tip,  but  also  the  effect  of  transition 

from  plastic  to  elastic  zones  in  V-V  enclosed  by  I\ 

e  234 

Now  we  shall  considi  r  the  physical  meaning  of  the  AT  integral  introduced 
in  Eq.  (93).  First,  we  1  ite  by  the  definition  of  the  potential  AU,  j 

1  T 

t  :Ae  +  AU  =»  t  :Ae  +  —  At  :  Ae 

=  (t+  ^  AtT) :  Ae  =  AW  (say).  (94) 

Thus  AW  is  the  total  stress-working  density  increment  during  the  time  interval 
t  to  (t+At).  As  shown  in [14]  and  [15]  we  may  write  AW  in  an  alternate  form, 
using  the  conjugate  variables  As  and  Ac,  where  As  is  the  incremental  second 
Piola-Kirehhof f  stress  such  that  (t+A*)  is  the  total  second  Piola-Kirchhof f 
stress  at  time  (t+At)  as  referred  to  the  conf iguration  at  t.  Thus: 
l  T 

AW  =  (t+  j  At  )  :  Ae 

l  T 

=  <t+  2  Af)  :  Ac  +  t  :  (A<-  .Ae) 
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=  (t+  2  As)  :  Ae.  (95) 

Recall  chat  t  and  As  are  both  symmetric,  and  the  symmetric  Ae  is  the  incremental 
strain  between  times  t  and  t  +  At  as  referred  to  the  configuration  at  t.  Using 
the  additive  decomposition: 

Ae  =  Ae  +  Ac  (96) 

-e  -p 

where  the  subscripts  e  and  p  refer  to  "elastic"  and  "plastic"  parts  respectively, 
we  see  that: 

t  :  e  +  AU  =  AW  =  (t+  ■—  As)  :  (Ae  +Ae  ) 

-  -  l  ~e  ~p 

=  AW  +  AW  .  (97) 

e  p 

Where  AW  is  "elastic"  part  of  incremental  stress-working  density  (per  unit 
e 

volume  at  time  t)  and  AW  is  the  plastically  dissipated  portion  of  incremental 

P 

stress-working  density.  Let  AK  be  the  incremental  kinetic  energy,  between  t 

and  t  +  dt,  per  unit  volume  in  b .  It  is  seen  that: 

AK  =  o  a  .  Au.  (98) 

t-t  - 

To  simplify  matters,  let  us  consider:  (i)  a  "mode  I"  2-dimensional  problem, 

wherein  the  crack-faces  are  free  from  ..ny  applied  i r act ions  or  other  constraints, 

(ii)  the  body  forces  to  be  zero,  and  (iii)  the  crack  is  along  the  x^  axis. 

Consider  a  volume  in  this  2-dimensional  problem,  which  is  enclosed  by  the 

contour  f  as  in  Fig.  2  .  Let  AY  ,  AW*  znd  AK*  be  the  incremental  potential 
2J4  111 

of  external  loads,  incremental  internal  energy,  and  incremental  kinetic 
energy,  respectively,  of  the  volume  V*.  It  is  seen  that: 


/*  =  f  AWdv  =  f 

1  Jv*  Jv 


[ t  :  Ae  +  AU  ]  dv 


:*  =  f  AKdv  =  /"pa 

1  Jv*  Jv  ^ 


(100) 
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4  JS 


Au. 


and 


Aip 


i  ■-  fr  l/V'ji 


+  At 


ji 


(y 


i)]  dpi| 


ds 


(101) 


234 


X 


Au 


r  +r  r 

‘234  12+  54 


(7nj  [Tji +  Atji  v1  dui! ds 


since  the  crack  faces  are  traction  free.  To  simplify  the  discussion  further. 


we  shall  consider  only  the  case  of  a  stationary  crack,  of  length  "c^".  The 
combined  incremental  potential  energy  between  t  and  t+dt  for  this  cracked  body  is 


AEX  =  -A<|i  -  AW*  -  AK* 


(102) 


Now,  let  us  consider  a  second  cracked-body  identical  in  geometry  to  the 
one  above,  but  with  a  crack  of  length  "c^+dc^".  Let  the  second  cracked  body 
be  subjected  to  a  load  history  that  is  identical  to  the  one  to  which  the  first 
body  has  been  subjected  to.  In  the  first  body,  let  the  cartesian  coordinates 


measured  from  the  crack-tip  be  *  Xl”*ci’  and  i’2*X2'  We  consider  in  the 

,2  • 


second-body,  a  contour  ^  which,  with  reference  to  the  boundaries  of  the  body. 


is  identical  to  the  contour  T it-  the  first  body,  in  other  words  one  may  say. 


equivalently,  that  we  ire  considering  identical  domains  V*  in  both  the  bodies 


except  that  the  leng^is  of  cracks  enclosed  by  the  contours  and  ^234’ 


respectively,  differ  by  an  amount  dc^.  Let  be  the  cartesian  coordinates 
centered  at  the  crack-tip  in  the  second  body:  =  c;^-dc^;  ^2  ~  ^2  *  X2’  Since 


both  the  bodies,  differing  In  crack-lengths  by  dc^,  have  been  subject  to  iden¬ 


tical  load  histories  until  time  t,  one  may  assume: 

=  T.  ,(r  )  ;  a_  =  a  (C,  ) 


in  body  1: 


ij 


cijvV  •  it  -  VV  ) 


Au  =  Au(r,^)  ;  AW  =  AW((;k) 


Tij  =  ^(rk);  ^  =  ^t(rk} 


in  bodv  2: 


Au 


Au(a')  ;  AW  =  AW  ( r, '  ) 
“  k  k 


(103) 
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The  combined  incremental  potential  energy  in  body  2  would  be: 


AE2  =  -A*  -  AW*  -  AK* 


(104) 


where , 


AW*  =  /  AW(c')  dv 

1  Jy*  k 


(105) 


AT2  =  Jv*  °t^t  aky  *  (Ck> 


(106) 


^i(5k> 


Aii»2  = 


Fr  +r  +p  » J 

234  12  54  0 


If  nj  [Tji(t;k)+  Atji(vii)1  dwijdS  (107) 


Integrals  in  Eqs.  (105-107)  are  identical  to  those  in  Eqs.  (99-101),  respectively, 
except  that  the  limits  of  integration  which  are  w.r.t.  5'  in  Eqs.  (105-107)  are 
different  from  those  w.r.t.  r, ^  in  Eqs.  (99-101).  By  using  arguments  analogous 
to  those  employed  earlier  in  connection  with  finite  elasticity,  it  can  be  shown 
that: 


AW*  -  AW*  = 


“/r  " 


^dCj ' Wds. 


(108) 


Note  that  singularities,  if  any,  in  AW  near  the  crack-tip  have  been  permitted 
in  both  bodies.  Since  the  crack  is  assumed  to  be  stationary,  one  may  assume 
that  neither  a  nor  Au  are  singular,  and  write: 


AT*  -  AT*  = 


*  =  -  r  ~L 
1  A*  3r’i 


(pa^ . Au) dc^dv 


(109) 


Au(c^)  =  Au  [  (r^-dc^)  ,r>2]  =  Au(r,R)  -  dc^ 


(110) 


Finally, 


A -  Aip 


p  3Au. 

i-+Jr  +r  +r  ”j  +  (TJi+ltJt>  ldeids-<lu> 

234  12  45  X  1 


Thus,  the  difference  in  incremental  total  energy,  in  the  time  interval  t  to 
t+dt,  between  the  two  bodies  is: 
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(112) 


AE2  -  AE 


i  "Jr 


, 


3Au_, 


n.  [ 


dc. 


234"+ri2+r45 


3  '34.  i 


Au,  +  (x.j+At..) 


/  " 


AWdv  + 


234 


ji  ji 

/  (pa.Au)dv. 

Jy*  d^i 


■]  ds 


However , 


f  n  ^  Au.  *  f 

J r^,  +  v^+r,z  j  ci  Jv 


Au, ) , .dv 


234  12  ‘45 


V* 


j i ,1  i  * j 


/  t (T ^ n-  ,)»,  Au.+t.  .  4n.  ]dv  (a13) 

Jv*  Ji,3  1  i  jx,loui,j 


In  writing  Eq.  (113),  it  is  assumed  that  is  non-singular  near  the  blunted 


crack-tip  in  an  elastic-plastic  material,  for  which  stress  saturates  to  a 
finite  value  even  for  large  strains.  It  is  recognized,  that  at  time  t  the  linear 
momentum  balance  condition  is: 


T  •  ,  j  p  (s_.  ) 


i  t 


‘  ji,  j 

or  (Tji^)^  =  MV,1’!* 


(114) 

(114a) 


ds 


dc. 


Using  Eqs .  (114a)  and  (113)  in  Fq.  (112),  it  is  seen  that: 

3Au, 

t^AW  -  n  (t^+At  ) 

'r, 

:  J4 

[~T,{  ,AU,  .  +  p(a4)t. 

ty* 

Comparing  Eq.  (115)  with  the  extreme  left  hand  side  of  Eq.  (93),  it  is  seen 
that 


ae2  '  AEi =  /  tniAW  "  nj(Tji+Atji)  -^7] 

r234 

r  3Aut 

+  L  [-TJi,laul,3  +»(“i>t— i 1  dv 


(115) 


ae2  ae^  ■  Al^ 


dc. 


where  Al^  is  the  x  component  of  the  path-integral  vector  defined  in  Eq.  (93).  Thus,  j 
in  general,  it  appears  that  the  AT  integral  of  Eq .  (93)  has  the  meaning:  it  is  the 
difference  in  combined  incremental  potential  energy  (ner  unit  crack  length  difference) 


in  the  time  interval  t  to  t+dt  of  two  bodies  which  are  identical  in  shape  and  load 
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history,  but  differ  in  crack  lengths  by  'dc^'. 

Note  that  in  the  derivation  of  Eq.  (115),  no  assumption  has  been  made 
concerning  any  loading/unloading  conditions  within  V*.  We  now  consider  certain 
cases  of  cracks  in  rate-sensitive  inelastic  materials. 

Rate-Sensitive  Inelasticity: 

A  rate-sensitive  constitutive  law  of  considerable  generality,  as  given 
by  Perzyna  [11], can  be  written  for  finite-strains,  when  an  associative  flow-rule 
is  used,  as: 


5a  =  *<♦<«>!£ 


(116) 


where  <>  denotes  a  specific  function,  such  that  <4>  =  4(F)  for  F>0,  and 
4=0  for  F<0.  The  parameter  y  is  called  the  fluidity  parameter  and  £a  is  the  in¬ 
elastic  strain  rate.  Various  forms  o'  4  were  reviewed  in  [11],  For  the  Hencky- 
Mises-Huber  type  yield  criterion,  one  can  define  V  to  be: 


F  =  [3J„(a')P  -  F  o  -  r  =  (~c':o'r  ~  F 
2  -  o  eq  o  2~  o 


(117) 


where  a  is  the  equivalent  Kirchhoff  stress.  A  simple  choice  for  4(F)  can  be: 


><F)  =  F  . 


(118) 


If  in  Eq.  (117),  one  sets  F^O,  the  viscoplastic  strain-rate  relation  (116) 
leads  to: 


.a  3  ,  ,n-l  , 
t,.  =  rr(o  )  o', 

ij  2'  eq  ij 

Defining  the  "equivalent"  creep  strain  rate  as: 

.  ,2  .a  .a  . V 

e  =  (r  £..e.  J 
eq  3  i]  ij 

it  is  seen  from  Eq.  (119)  that 


(119) 


(120) 
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eq 


*<°eq> 


which  is  the  well-known  power-law  for  steady-state  creep.  If  d  is  the  co- 
rotational  stress-rate  (=Aa/At)  the  constitutive  law  can,  in  general,  be  wri- 
ten  as: 


5  =  L  :(£-ea)  (122) 

-  ;e  -  -  , 

where  L  is  the  tensor  of  instantaneous  elastic  moduli!,  and  ia  is  the  strain- 
~e 

rate  under  viscoplasticity,  or  creep,  respectively,  as  given  in  Eq.  (116)  or 
(119).  Thus, 

t  +At 

Ao  =  Le:Ae  -  Le:  f  jadt.  (123) 


From  Eq.  (123),  it  is  seen  that  a  potential  aV  for  an  elasto-viscoplastic 
solid,  analogous  to  that  in  Eq.  (73),  can  be  constructed  [  16  ] .  Using  Eq. 

(77),  a  potential  AU  for  At,  in  the  case  of  viscoplasticity  or  creep,  analogous 
to  that  in  Eq.  (79),  can  be  constructed.  The  main  points  to  be  noted  however, 
are:  f i )  in  the  case  of  an  elasto-viscoplastic  solid  AU  will  be  an  explicit 
function  of  position  of  the  material  particle,  since  <4>  (F)>  is  non-zero  or 
zero  (depending  on  F>0  or  <0)  at  a  given  location;  (ii)  in  the  case  of  creep 
AU  will  not  be  an  explicit  function  of  location  of  the  material  particle. 

From  the  above  discussion,  it  is  evident  that  in  the  case  of  cracks  in 
elasto-viscoplastic  solids,  an  incremental  integral  vector,  AT,  entirely 
equivalent  to  Eq.  (93)  can  be  defined.  Moreover,  it  is  also  evident  that  the 
AT  integral  in  the  case  of  viscoplasticity  is  path- independent  in  the  same 
limited  sense  as  discussed  earlier  in  connection  with  rate-independent  elasto- 
plasticlty . 

On  the  other  hand,  in  the  case  of  creep,  described  by  a  constitutive 
law  of  type  (119),  since  AU  is  not  an  explicit  function  of  position  of  the 
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material  particle,  we  have  the  result,  analogous  to  the  case  of  finite  elastic¬ 
ity,  that  R=0  in  Eq.  (93). 

Thus,  in  the  case  of  creep,  we  have,  for  instance  in  a  2-dimensional  case, 
the  following  result  instead  of  Eq.  (93): 


L 


234 


[n(x. ’Ae+AU)  -  n  .  <(x+At)  .  Ae>]ds  +  Lt 
l  '  '  e+0 


/ 

JM-\ 


[-V  x:Ae 
— t  ~ 


-  p  (_f-a)  .  Ae]dv  + 


J *  n+  [(x:Ae 


Ae+AU)  -  (x:Ae+AU)  ]ds 


J *  [(_t.  Ae)  +  +  (t_.Ae 

ri2 

J"  [n(x:Ae 


12 

Ae)  ]ds 


/ 

12 


[<(x+At)  .  Ae>  -  <(x+At).Ae>  ]  ds 


=  Lt 

e-*- o  -'r 


Ae+AU)  -  n  .  < (x+At )  .  Ae>]ds 


=  AT^  (124) 

v  lere  ,  the  incremental  vector  for  the  case  of  a  crack  in  a  creeping  material, 

is  once  again,  bjf  definition,  the  non-zero  limit  (in  the  limit  e-*-0)  of  the 

< ontour  integral  on  F  as  in  Eq.  (124).  The  far  field  evaluation  of  AT 
°  e  -  — c 

should  be  {.  irformed  as  per  the  contour/ volume  integrals  appearing  on  the 

extreme  L.H.S.  of  Eq .  (124).  Once  again,  it  is  evident  that  since  AU  is  not 

an  explicit  function  of  location,  the  integral  appearing  on  the  extreme 

L.H.S.  of  Eq.  (124)  is  strictly  path-independent  (cventhough  the  enclosed 

volumes  V-V  change).  Further  it  is  seen  that  the  AT  integral  in  the  case 
e  — c 

of  creep,  characterizes  the  severity  of  crack-tip  conditions. 

Consider  the  application  of  F.q .  (124)  to  the  special  case:  (i)  body 

forces  including  inertia  are  negligible,  (ii)  the  crack-faces  are  free  from 

any  prescribed  conditions,  (iii)  the  crack  is  aligned  with  the  axis  and  that 

we  are  interested  only  in  the  x  component  of  the  vector  AT  .  Under  these 

1  ~ c 

restrictions,  Eq.  (124)  becomes: 
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X 


3Au, 


234 


[n^rrAe+AU)  -  — 


]ds  -  Lt 

1  e-0 


If  3Tu  3Aut 

(Jv_v  a*L  3X. 


dv! 


=  Lt 
e-0 


y*  [n^ (t  :Ae 


3  Au . 


:Ae+AU)  -  n.  (r..+At..)  — — 
3  Ji 


— ]ds 


5  (4T1>C- 


(125) 


Once  again,  it  can  be  shown  that  (AT  )  is  the  difference  in  incremental  total 

1  c 

energies  during  t  and  t+At  of  two  bodies  identical  in  shape  except  differing 
in  crack-lengths  by  'dc^*,  both  of  which  have  been  subject  to  identical  load 
histories  upto  time  t. 

We  observe  that  Landes  and  Begley  [10],  and  Goldman  and  Hutchinson  [9]  in¬ 
troduced  a  parameter  c*  defined  as: 

3u. 


r  du. 

c*  =  /  [w*dy  -  T± 


(126) 


to  correlate  creep  crack  growth  rates  in  "Mode  I"  where  r  is  a  far-field  contour 


[analogous  to  r ^  in  Eq.  (125)]  and 


W*  = 


tf  d£  .  . 


(127) 


ij  iJ 

we  note  that  the  above  parameter  was  introduced  in  [9,  10]  based  primarily 

on  the  observation  of  the  essential  similarity  between  the  pure  power-law 

hardening  type  (deformat ion-theory-of  plasticity)  constitutive  relations  between 

the  strain  and  stress  on  the  one  hand,  and  the  steady-state  creep  relation  of 

the  power-law  type  between  the  strain-rate  and  stress  on  the  other. 

The  present  parameter  (AT^)^,  of  Eq .  (125),  however,  accounts  for  elastic 

g  c 

as  well  as  creep  strains  (using  the  decomposition,  Ac  =  Ae  +Ae  where  super¬ 
scripts  e  and  c  stand  for  "elastic"  and  "creep"),  finite  deformations,  and  is 


valid  even  under  non-steady  creep  conditions.  The  path-independent  nature  of  (AT,) 

1  c 


is  in  Eq.  (125)  is  evident  from  previous  discussion.  With  the  definition  of 


(AT^c  as  the  limit  e-K)  of  the  contour  integral  T^,  the  local  nature  of 

(AT,)  ,  i.e.,  its  role  as  a  parameter  quantifying  the  crack-tip  conditions  is 
1  c 

evident. 


Under  steady-state  creep  conditions,  it  is  seen  that  (AT^)c  of  Eq.  (125) 

g 

reduces  to  (AT,)  where, 

1  c 


/•  3Au. 

J  lni(I:As)“  nj(Tji)  J77] 


(128) 


or  (Vc  =  et-O  /  tnl(I:-4)  -  nj(Tji)  ds 


(129) 


on  the  other  hand,  C*  of  Eq.  (126)  can  equivalently  be  written  as: 


r-  3u. 

J  [W*dy  -  Tt  -—]  ds 


(130) 


where 


f 


a.  ,d£  . 
ij.  iJ 


(130b) 


Thus,  e'en  under  steady-state  conditions,  the  present  (t, )  differs  from  C* 
of  [9,10],  in  that;  while  i:e  occurs  in  the  former,  W*  occurs  in  the  later. 

It  is  noted  that  while  t:e  has  the  physical  interpretation  as  the  stress-working 
rate,  W*  is  simply  a  mathematical  potential  for  o^_.  under  steady-state  creep 

•  g 

conditions.  Thus  while  (T, )  still  has  a  physical  interpretation  as  discussed 

I  c 

earlier,  C*  does  not,  in  general,  have  a  physical  meaning.  These  and  other 


topics  relevant  to  the  use  of  (AT^J^  and  C*  are  more  elaborately  discussed  in 
[17]. 
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